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Abstract 

We derive an approximate analytic formula for the ground-state energy of the charged anyon 
gas. Our approach is based on the harmonically confined two-dimensional (2D) Coulomb anyon 
gas and a regularization procedure for vanishing confinement. To take into account the fractional 
statistics and Coulomb interaction we introduce a function, which depends on both the statistics 
and density parameters (y and r s , respectively). We determine this function by fitting to the 
ground state energies of the classical electron crystal at very large r s (the 2D Wigner crystal), and 
to the Hartree-Fock (HF) energy of the spin-polarized 2D electron gas, and the dense 2D Coulomb 
Bose gas at very small r s . The latter is calculated by use of the Bogoliubov approximation. 
Applied to the boson system (v = 0) our results are very close to recent results from Monte 
Carlo (MC) calculations. For spin-polarized electron systems (v = 1) our comparison leads to a 
critical judgment concerning the density range, to which the HF approximation and MC simulations 
apply. In dependence on u, our analytic formula yields ground-state energies, which monotonously 
increase from the bosonic to the fermionic side if r s > 1. For r s < 1 it shows a nonmonotonous 
behavior indicating a breakdown of the assumed continuous transformation of bosons into fermions 
by variation of the parameter v. 



I. INTRODUCTION 

The jellium model consisting of a homogeneous electron (or fermion) gas charge- 
neutralized by a uniform positive background, is serving as a realistic model to describe 
generic properties of metallic systems, which extend in two (2D) or three (3D) spatial 
dimensions [2]. Especially, the exchange and correlation contributions to the ground state 
energy per particle in dependence on the particle density (quantified by the Wigner-Seitz 
radius r s ) and the degree of spin polarization are required in applications of the density 
functional theory (in its local density approximation) to inhomogeneous fermion systems 
(e.g., nuclei, molecules, solids) [3j]. Analytic formulas of these contributions for 2D and 3D 
systems are available from interpolations between known data from numerical or analytical 
calculations for special values of the system parameters Although frequently in 

use, such formulas are not free of shortcomings and may lead to misleading results Here, 
the special topology of 2D systems provides a unique opportunity to derive an approximate 
analytic formula for the ground state energy, which has not been exploited so far. 

n 

The 2D topology allows for fractional exchange statistics [8j , characterized by a continuous 
parameter u, that may attain values between (for bosons) and 1 (for fermions). Particles 
with < v < 1 are generically called anyons []l3| • The quasiparticle excitations in the 



fractional quantum Hall regime 
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13j and in certain quantum magnets [14) can be 



described using the anyon concept. 

2D electron systems are realized in semiconductor heterostructures and layered metallic 
systems such as cuprate superconductors. Their ground-state properties are subjects of 



fundamental studies (see Ref. 15( and references therein). To obtain accurate estimates 
of the ground-state energy, numerical simulations have been employed 4], [5], which are 
restricted, however, to special values of the particle density (expressed by the interparticle 
distance r s ) or the degree of spin polarization. An analytic expression in terms of r s exists 
for the 2D Wigner crystal (r s ^> 1) [16], which represents a classical system. All these 
investigations refer to 2D electron systems and are, as such, based on the fermion character 
of the particles. It is interesting to note that bosonic Coulomb systems have found so far 
only little attention 17|, which may be due to lacking realization. This seems to be even 
more the case for the anyon aspects of the 2D Coulomb gases, while for the 3D Coulomb 
Bose gas a closed-form expression for the ground-state energy as a function of r s was derived 
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by use of the Bogoliubov approximation [3]. The aim of this paper is, by making use of the 
anyon concept, to establish a link between the analytic results, which exist for the ground- 
state energies of the 2D Coulomb Bose gas (derived here in close analogy to Ref. [181]). of 
the Wigner crystal [16], and of the high-density 2D spin-polarized electron gas. 

Previously, we have derived an approximate analytic expression for the ground state 
energy of N charged anyons confined in a 2D harmonic potential [3]. This was achieved 
by using the bosonic representation of anyons and a gauge vector potential to account for 
the fractional statistics, which allowed us to work with the product ansatz for the TV-body 
wave function. A variational principle has been applied by constructing this wave function 
from single-particle Gaussians of variable sha pe. As in many other perturbative treatments 
of anyons in an oscillator potential (see Ref. [l9| and references therein) our expression for 
the ground state energy had a logarithmic divergence connected with a cutoff parameter for 
the interparticle distance. Making use of the physical argument (see Ref. [l2j) that for v ^ 
this distance has to have some finite value, we have regularized the formula obtained for 
the ground state energy by an appropriate procedure that takes into account the numerical 
results for electrons in quantum dots in the case with Coulomb interaction. 

The ground state of anyons in a harmonic potential has been studied by Chitra and 
Sen [20y , especially for the limit of an infinite number of particles, yet without Coulomb 
interaction. In the present paper, we make use of previous work for the harmonically confined 
2D anyons {19I . I20I ] to derive an approximate analytic expression for the ground-state energy 
of the homogeneous 2D anyon gas. This is done (for all values of the statistic parameter v) 
by flattening out the confining potential with a simultaneous increase of the particle number 
N, but fixed areal density, to obtain the infinite size system, i.e., the thermodynamic limit. 
It is achieved by redefining the strength u>o of the harmonic potential such that it vanishes 
with increasing N as 1/iV 1 / 2 . With respect to the dependence on the statistic parameter 
we have to distinguish between the cases without and with Coulomb interaction. In the 
former case, we introduce a function of z/, which is fitted to known analytic expressions 
at the fermion or boson end. For the case with Coulomb interaction, one has to warrant 
also charge neutrality. This is achieved by introducing the interaction with a positively 
charged background disk, which slightly modifies the interaction term in the formula for the 
ground-state energy. For this CclS6 ; clS the dependence on v may now be mixed with that on 
the density parameter r s , we assume a function, which depends on both parameters. This 
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function is fitted to the analytic expressions for the ground state energies of the classical 
Wigner crystal (for large r s and independent of v) and of the 2D Coulomb Bose gas (for 
small r s and v = 0). For the fermion case [y = 1), we fit this function to the analytic 
expression for the HF energy, which is a high-density limit (r s — > 0). 

In Sec. [Ill we introduce the principal approach for deriving an analytic formula for the 
ground state energy and demonstrate the regularization to obtain the thermodynamic limit. 
In Sec. II III we derive an analytic expression for the ground-state energy of the 2D Coulomb 
Bose gas using the Bogoliubov approximation. In Sees [TV] and [V] the derivation of the 
approximate analytic formula for the ground state energy of the 2D Coulomb anyon gas is 
completed and the results will be discussed by comparing with numerical results from the 
literature. We summarize and conclude with Sec. IVII 



II. PRINCIPAL APPROACH AND HARMONIC POTENTIAL REGULARIZA- 
TION 



Let us start by considering N noninteracting spinless anyons of mass M confined in a 2D 
parabolic potential described by the Hamiltonian 



I N r 2 



(1) 



Here r*j and Pi are the position and momentum, respectively, of the zth anyon in 2D and 



N 



e 2 x n 



(2) 



is the anyon gauge vector potential 
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22j with fij = r*j 



fj , and e z the unit vector normal 



to the 2D plane. The parameter v accounts for the fractional statistics of the anyon: it 
varies between v = for bosons and v = 1 for fermions. When later on considering the 
anyons as charged particles, we add their mutual interaction to define the Coulomb anyon 
problem. 

In Ref. we have outlined a variational procedure for the ground-state energy of 

interacting anyons confined in an oscillator potential with characteristic energy Hujq. Starting 
from the bosonic end [y = 0) we achieved, after regularization of a logarithmic expression 
by means of a cutoff parameter for the particle-particle interaction, approximate analytic 
formulas in terms of iV and v. For the noninteracting anyon system we found 

E (N, v) = hu NAf 1/2 , (3) 



while for the interacting anyon system it is given by 



E (N, v) 



M + X 2 



2M 



with 



and 



X = (A + B) 1/2 + {-(A + B) + 2{A 2 -AB + B 2 ) 1/2 ] 1/2 



A = [M 2 /128 + ((AT/12) 3 + (M 2 /128) 2 ) 1/2 ' 
B = [M 2 /128 - ((AT/12) 3 + (M 2 /128) 2 ) 1/2 

In these expressions we use Af = 1 + v(N — 1) and 

/71V/2 N - 1 L 

M 



1/3 
1/3 



(4) 



(5) 



(6) 



(7) 



.2) 2 a B 

where L = (h/Muo) 1 ^ 2 is the oscillator unit length, and as the Bohr radius. 

In order to obtain the corresponding expressions for the homogeneous 2D anyon gas, we 
flatten out the parabolic confining potential while increasing the number N of anyons, but 
keeping the density p = N/S = l/nr 2 constant, i.e. we perform the thermodynamic limit 
while making the confining potential disappear. Here nr 2 . is the area of the jellium disk 
carrying the positive countercharge corresponding to the mean particle distance ro = asr s 
expressed in units of the Bohr radius by the dimensionless density parameter r s . 

Without Coulomb interaction and in the case of fermions {y — 1), the ground state energy 
of the homogeneous 2D electron system of density p is determined by the Pauli exclusion 
principle. It is given by 

E (p) = irh 2 pN/M, (8) 
while from Eq. ([3]) we have (see also Ref. [2o| for v = 1 and iV — > 00) 

) = huo N 3 / 2 . 



E (N,v 



(9) 



In the thermodynamic limit both expressions have to become identical and we obtain the 
relation 

uj (N) = 7rhp/(MN 1/2 ) , (10) 

which means that, in fact, the thermodynamic limit (N —>■ 00) is obtained for the vanishing 
parabolic confining potential. We extend this consideration for the fermionic limit, v = 1, 
to the general anyon case, v ^ 1, by assuming instead of Eq. (jSJ) the relation 

E Q {p,v) =7rh 2 pN<p(u)/M, (11) 



where the function 4>{v) is still to be determined under the constraint <f>(u = !) = !. This 



form is motivated by the fact that close to the bosonic 
of the infinite anyon gas depends linearly on v 23l . 



imit {v ~ 0) the ground-state energy 
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25j. We replace this particular 



dependence by <fi(v). Consequently, we have to change Eq. ([TO]) into 

uj (N,u) = nhpf(u)/(MN 1 / 2 ) (12) 

with another unknown function f{y) and the constraint f(y = 1) = 1. Thus we assume the 
vanishing of Uq according to 1/iV 1 / 2 also in the general anyon case v ^ 1. As it turns out 
<p{v) is determined by f(v). For the free (noninteracting) anyon gas both unknown functions 
would be the same. However, in our generalization to the charged anyon gas this will not 
be the case anymore. 

In the thermodynamic limit, and including the Coulomb interaction, the parabolic con- 
finement (caused by geometry and electrostatics of the compensating charges) has to be 
replaced by the jellium contribution, which for a disk of radius Rq (containing iV counter- 



charges) gives a potential energy contribution [26] 



V(fk) = -p[ ■ (13) 

Js \r k - r\ 

Here S = ttRq is the area of the jellium disk for N charges and we have Ro = N 1 l 2 r$. 
For v = 1, we may now identify the characteristic length L of the oscillator with the mean 
particle distance to obtain L = {%/ 'Mojq) 1 ^ 2 = iV 1//4 r , and find for N ^> 1 the relation 
r < I <C The generalization to the anyon case, v ^ 1, is possible by dividing those 
lengths by / 1/2 (z/). 

The strength of the particle interaction is characterized by the density parameter r s (in 
2D the radius of the Seitz circle), and we may generalize our analytic expression ( TTTT) for 
the ground state-energy to the form 

Eo(v, r s ) = nh 2 P N(P(v, r s )/M . (14) 

Thus we replace the unknown function (j){y) by (j)(u,r s ), which now depends on the two 
system parameters v and r s and takes into account the Coulomb interaction and the statis- 
tics. In the high density limit (r s — ► 0) it becomes the function 4>(u) introduced for the 
noninteracting system. The nontrivial case of r s — > and v — > will be considered in Sec. 

El 
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In a similar way, we have to generalize f(u) to f(u,r s ). Now these functions depend in a 
complicated way on each other. Thus we have established the extension from the confined 
to the homogeneous 2D anyon system. In Sec. |V]we determine the function f(u,r s ). 



III. 2D COULOMB BOSE GAS AT HIGH DENSITIES 



The 3D Coulomb Bose gas problem has been treated by Foldy [18j in the high density 
limit by applying the Bogoliubov approximation [27]]. Originally, Bogoliubov considered a 
ow-density system of bosons interacting with short-range forces, but it was shown in Ref . 
l^l that the Bogoliubov approximation can also be reliably used for bosons with long-range 
Coulomb interactions in the high density limit. In this section we use the same strategy for 
the 2D Coulomb Bose gas to obtain an analytic expression for the ground-state energy. 

According to Ref. 18J, the Bogoliubov approximation is valid if almost all particles are 
in the zero momentum state, i.e., (N — N Q )/N tends to zero when r s — > 0, where Nq is the 



181 ]. for the 2D Coulomb Bose 



number of particles with zero momentum. Following Ref. 
gas, this ratio takes the form 

N 4vriV Jo \ E k J K J 

where £ k = % 2 k 2 / (2M), V k = 2ire 2 / (Sk) is the Fourier transform of the 1/r potential in 2D, 
and 

E k = [(S k + NV k ) 2 -N 2 V 2 ] 1/2 
is the dispersion of collective excitations. By introducing the variable 

h 2 k 3 s \ 1/6 



(16) 



^ {^nMNe 2 J 



a 3 B r 2 k 3 \ ^ 6 



(17) 



Eq. ( fToT) takes the form 



(N - N ) r 2 J 3 



N 



2V3 Jo 



e 6 + 2 
(e + 4)V2 



After evaluation of the integral one obtains 

(n-Nq) r(|)r^ 7 



eZ r 2/3 



0.701091r? /3 



(18) 



(19) 



N 2^ 

which tends to zero for r s — > 0, thus showing the validity of the Bogoliubov approximation 
also for the 2D Coulomb Bose gas. 



The 2D analog of the ground-state energy given in Ref. 18J is 

Q roo 

E — — (E k -£ k - NV k )k dk , (20) 

47T JO 

and can also be written in terms of the variable £ (in Ry=e 2 /(2ob) units) 

I ie(e+i) i/2 -e-md^ (21) 

which after evaluation of the integral yields 



E. - 01/3--2/3 



| = -CB G r 7 2 / 3 , (22) 



where cbg = — 3 ^ 6 = 1.29355. The same result has been obtained also in Ref. [28[ by 
applying the hypernetted chain approximation. 

Thus we have an exact analytic expression for the ground-state energy per particle for 
the 2D Coulomb Bose gas valid at high densities. It will be used in Sec. |V] together with 
the known expression for the 2D Coulomb gas in the low density limit (the 2D Wigner 
crystal) to derive a form for the unknown function f(v,r s ) introduced in Sec. [TT1 We note 
in passing that the 2D (classical) Wigner crystal is independent of particle statistics, thus 
its properties do not depend on v. 

It is interesting to note that the spectrum of collective excitations of the 2D Coulomb 
Bose gas, Eq. ffl"6l) . can be cast in the form 

E k = [h 2 27ie 2 pk/M + (h 2 k 2 /(2M)) 2 } 1/2 , (23) 



which for small k represents the 2D plasmon dispersion 



29J and approaches for large k the 



free particle dispersion. In contrast to the 3D system this spectrum has no gap at k — 0. 



IV. ANALYTIC EXPRESSION FOR THE GROUND STATE-ENERGY IN THE 
THERMODYNAMIC LIMIT 

With the considerations of the previous sections we may now derive the wanted analytic 
expression for the ground-state energy of the charged anyon gas. The system Hamiltonian 
is given by 

N 

1=1 




+ V(?i 



(24) 
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where (for iV — > oo) V(fi), the interaction energy of the ith electron with the jellium 
background, plays the role of the confining potential considered in Ref. jl9j ]. 

From our calculation for parabolically confined interacting anyons we know the contribu- 
tions of the individual terms of Eq. (1211) to the expectation value of the ground state energy 
except for the potential energy V{fi). As will be demonstrated in this section, V{fi) modi- 
fies the expectation value originating from the Coulomb interaction in the case of A" — > oo. 
Namely, we will show that, differing from Eq. ([7]), A4 will now be proportional to N 1 ^ 2 . 
This allows us to obtain in Sec. |V]the thermodynamic limit of the ground-state energy, Eq. 
(BJ, as a complicated function of the system parameters v and r s introduced by M. and its 
dependence on /(z/, r s ). 

Before this is done, we have to calculate the mentioned modification of the contribution 
of the Coulomb interaction to the ground-state energy by V(fi). For this we adopt the 



approach of Fisher et al. [30j who used the Ewald method. They started dividing the plane 



into square cells, each cell containing the same number of pointlike charges (electrons) and 



the corresponding 



ellium background. We use their expression for the interaction potential 



. 3 



(Eq. (A2) of Ref. 

2tt _ erfc(SY)e^ e rfc(|f - K\/2t) Air^H 

$(r) = Z^ s + £T- W 

Here, TZ = C(l x ,l y ), and S = 2n{l xi l y ) / C (with l x ,l y = 0,±1,±2,...) are the direct and 
reciprocal lattice vectors, respectively, with lattice constant £; erfc(x) = 1 — erf(x) is the 
complementary error function, and t the Ewald parameter for optimal convergence of the 
lattice sums. 

The interaction energy of the 2D Coulomb gas is given by 

,2 



1 N 



V(f u ...f N ) = 



Hn - fj) - 5 k 



J ) l 3 i -» -» 

i,j=i i i'» 'ju 



(26) 



where the self-interaction is explicitly subtracted. Using the Gaussian variational wave 
function of Ref. Jjj] to compute the expectation value of the interaction energy, we need to 
evaluate 

]T J d\ d 2 rj e- ar ^ e- ar 1$(f i - f 3 ) . (27) 

Consider now the individual contributions to $(rj — fj) in Eq. ( 1251) . Typically, the Ewald 
parameter is taken to be t ~ C. For simplicity we assume t = C Then St = 2u:{l 2 x + l 2 ) 1 ^ 2 , 



— * — * 

\r — lZ\ — C\r/C — (l x , ly)\, and \f— 1Z\/[2t) = \r/C — (l x , l y )\/2 and we rewrite the expression 
for $(r), Eq. (125|) . in the form 



where 



A{r) 



E 



_ . erfc(r/2£) 1 

$ r = — - + -i r 



erfc(27r(e + ?) 1/2 )c^ , erfc(|f7£-(^g|/2)' 



(28) 



(/2+/2)V2 



|f/£ - (Za;,Zy)| 



47T 1 / 2 . 



(29) 



The length of the vector r*is in the interval < r < C We identify £ = L = iV 1//4 r , 
which implies r$ < r < iV 1//4 ro, and estimate the contributions to $(r) in the limits r <^ L 
and r — > L. For r/L <C 1 one has e* S r ~ 1 and one can neglect f/C in the second term of 
A(r) to obtain the leading terms of an expansion $(r*) w l/r + A/L. When r —+ L and with 
S -f— 2ii(ll + ly) 1 / 2 cos /3 , we find that the first term of A{r) depends only on l x , l y and the 
angle between 5* and r. Likewise, it is seen that the second term of A{r) neither depends 
on r nor on L. Hence we conclude with the relation 



$(r — > tq) 



iV 1//4 , (30) 
oo) the small r contributions to the inter- 



$(r -> L) 

showing that in the thermodynamic limit (N 
action energy become essential. 

Thus, assuming $(r) w 1/r, and taking into account that the original cell contains 
irL 2 / {Anr 2 ) = N 1 ^ 2 /4 particles, we obtain (after integration) for the expectation value of 
the interaction energy per particle 



N 



(31) 



Rewriting the result for E^ in h 2 /(ML 2 ) units, substituting into the expression for the 
expectation value of the total energy, and minimizing with respect to a as in Ref. [3] yields 
the expression for the ground-state energy, Eq. with Hlu = h 2 /(ML 2 ) and 

aN l / 2 L 



M 



a B 



(32) 



Following Ref. [19J we assume Af = 1 + v(N — 1). The constant a is fixed by fitting our 
ground-state energy to the exact numerical value of the 2D Wigner crystal ground state 
energy and we find a = —Cwc- 
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V. GROUND-STATE ENERGY OF 2D COULOMB ANYON GAS 



For sufficiently large N and by taking into account that L = iV 1 / 4 ?-^, we may use Ai = 
— cwcN 3 ^r s / f 1 ^ 2 ^, r s ) and J\f = uN (in this section the boson limit, v = 0, is understood 
as the limit v — > under the condition v ^> 1/N ) to write, starting from Eq. the 
expression for the ground-state energy per particle (in Ry units) in the form 



2K\ 



+ 



K 2 



K 



K 



x 



Here 



and 



K x = (K A + Kb) 1 ' 2 + [-(K A + K B ) + 2{K\ - K A K B + K 2 B f 



/211/2 



(33) 



(34) 



K A 
K B 



(35) 



K 2 /128 + {{v/l2f + (K 2 /128) 2 ) 1/2 ] 1/3 , 
K 2 /128 - {{v/l2f + (K 2 /128) 2 ) 1/2 ] 1/3 , 

with K = cwcTsj f x ' 2 (y-, ?"s), are obtained from Eqs. (jSJ) and (jSJ) for X and A and S, 
respectively, after substitution of M. and A/" by the expressions given above. Note that 
we have scaled the lengths by / 1 / 2 (z/, r s ) (see Sec. [II]) and replaced c^o according to Eq. 
(112p . Finally, we determine the function f(i/,r s ) by fitting Eq. ( 1331) to the known analytic 
expressions for the ground state energy of spin-polarized electrons and of the 2D Coulomb 
Bose g derived in Sec. IHII 

For the Bose gas [y = 0) we obtain from Eq. (|33|) 



£0(0,7%, 



2/3 



/ 2/3 (0, 



4/3 



(36) 



and find with Eq. fl22l for small r s /(0, r s ) ~ CBG r s/ c wc- F° r large r s , the ground-state 
energy does not depend on statistics and equals the energy of the classical 2D Wigner crystal 
Q, E wc = -2.2122/r a . This matches with Eq. (PJ if at low densities f(0,r a ) « r 1 / 2 with 

2 /3 

c iyc = 2.2122. For arbitrary r s , we connect these asymptotics by 



/(0,r s ) 



cIc^/Wc . 0.2r 2 Ln(r, 



+ 



(37) 



1,3/2 1/2/ 1 1 1 r 2 

1 + c BG rs l C WC l + r s 

For z/ 7^ and small r s the asymptotics of the ground-state energy, Eq. (I33p . has the form 

2/(f,r a ) f 1/2 civc^ 



0) 



7/ 



z/ l/4 / l/2 (z/;rs) 



(38) 
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the first term of which (for f(u,r s ) = 1) corresponds to the thermodynamic limit of the 



ground-state energy of anyons in 2D harmonic potential without Coulomb interaction 



201 ] . For large r s we obtain 



S (u,r s -> oo) = wcJ 1 "M -1 + ' K 4 y J . (39) 

The function f(i/,r s ) has to fulfill the following constraints: /(z/ = l,r s — 0) = 1 for the 
dense ideal Fermi gas, f(v,r s = 0) = z/ 1 / 2 for the ideal anyon gas close to the bosonic limit, 
and /(0, r s ) given by Eq. (|37[) for the 2D Coulomb Bose gas. The interpolating functional 
form 

SM « ^(r.).* + - fffi° + ^MfiM (40) 

1 + ci (r a )c B ' G 7V / Ch/c 1 + r s 

with c (r s ) = 1 + 6.9943r s + 22.4717r 2 and Ci(r s ) = 1 — e~ Va satisfies these constraints and, in 
addition, yields in the fermion case {y = 1) for the ground-state energy per particle the HF 



result 



31] . Moreover, for intermediate r s the logarithmic term in Eq. fj40l gives ground-state 



energies of the Bose gas lower than those of the spin-polarized electron system, as indicated 



in 



In Fig. 1 we show results for the ground-state energy per particle on the large scale 
1.0 < r s < 15.0. The upper four curves refer in descending order to the fermion case [y = 1): 



HF energies for spin-polarized electrons from Ref. 3jJ (open triangles), interpolated by Pade 



approximant MC data from Ref. {5] (on the given scale identical with those of Ref. jg]) 
for spin-polarized electrons (crosses), and our results from Eq. (133]) (on the given scale 
identical with those of Eq. ([39]) ) for spin-polarized electrons (closed circles). The lower 



two curves are for charged bosons (v = 0) and result from MC calculations of Ref. 17] 



(closed triangles) and from our Eq. (1551) (open squares). By star symbols we indicated the 
MC data 4j (without interpolation) obtained for some particular r s values. In Fig. 2 the 
corresponding results are depicted for 0.2 < r s < 1.5 and a larger energy scale. Here the 
difference in the data from Eqs. (!33l) and (!39|) (plus signs and closed circles, respectively) 
is clearly resolved for the smallest r s values. 

Let us first discuss the curves for the interacting Bose gas. As can be seen from Figs. 1 
and 2, our analytic formula (13"6T) for the ground-state energy yields results in close agreement 
with the MC data of Ref. [171], which are obtained numerically with much effort and represent 
upper bounds to the ground state. In this respect it is noteworthy that for r s > 2.0 the 



energies of Eq. (1361) are lower than those of Ref. [171 ] . 
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For the spin-polarized electron system we recognize, that the results of our analytic 
formula (Eq. fl33|) ) are lower than the HF data for all r s with larger deviations for small 
r s , while the MC data (which are close to the HF results) correctly approach the Wigner 
crystal limit (as our results). The curves obtained from our low-density expansion, Eq. 
( 13"9"j) . are very close to those of Eq. (I3"3"|) except for the lowest r s values, thus indicating to 
breakdown of this approximation (Fig. 2). The deviation from the HF and MC data for 
smaller r s can be made explicit by looking at the minima of these curves. By minimizing Eq. 
( 139]) with respect to f 2 / 3 {v,r s )/r^ 3 we obtain £ , min = — 3c^ c /(28z/) for v ^ 0. For the 
electron system [y = 1) this gives with cwc = 3.2903 the energy minimum £o,mm = —1.159. 
With f(u = l,r s ) from Eq. (I40p its position at r s m 0.7 is in close agreement with the 
result obtained from Eq. ([3"3"1) (as shown in Fig. 2). In contrast, the HF energy [31] for the 
spin-polarized electron gas 

Ehf=-,-^~ (41) 



r? 37rr 



s 



takes its minimum value Ehf ,min = —0.360 for r s m 2.36 compared with MC data of 
^MC,min = -0.393 at r s « 2.3. 

Here we have to consider that the validity of the HF approximation is limited by the 
demand that the leading (kinetic energy) is larger than the second term (exchange contri- 
bution). For the spin-polarized electron system this is the case for r s < 1.2, i.e., the HF 
minimum energy per particle is achieved in a density range, for which the HF approxima- 
tion does not apply. Likewise, we can discard the MC results for this intermediate density 
range because they are obtained with a method, which conceptually is close to the HF 
approximation. 

Results for the v dependence of the ground-state energy per particle, calculated for various 
fixed values of r s using Eq. ( |33|) . are presented in Fig. 3. It shows for r s > 1 a monotonous 
increase from the boson [v = 0) to the fermion [y = 1) end. This follows analytically also by 
taking the derivative of Eq. (1391) . valid for this range of r s . In contrast, for r s < 1, a minimum 
appears for an intermediate value uq, thus the boson energies are not the lowest ones. This 
minimum follows (for v < r s ) with f(v,r s — > 0) ~ v 1 l 2 b\ with h\ = 1 + 2.441472r s from 
Eq. (HOT) using the approximate Eq. fl3"9~|) . It is given by S 0imin = —(A/fybl^cwcbl^/rs with 
62 = 3/35, and occurs at z/ = b\J 4 cwc r s/b l i 2 . For r s = 0.2 this gives £o,?™n ~ —8.689 and 
u 0.08, while from Eq. (133!) we have £o,mm = —8.847 and u = 0.08. From this analysis 
we see that for r s — > the minimal energy decreases according to £o,mm ~ — l/ r s much 
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FIG. 1: Ground-state energies per particle vs density parameter r s for range 1.0 < r s < 15.0 from 
top to bottom: for fermions (y = 1) HF approximation (Eq. (|4ip and Ref. [3l|, open triangles), 
MC interpolation data (from Ref. [5j, crosses), and present results from Eq. ([39]) (closed circles) 
and Eq. (|53"j) (plus signs), and for bosons {y = 0) present results from Eq. (f5U|) (open squares) 
and MC data from Ref. [17] (closed triangles). MC data of Ref. [4] for some particular values of 
r s are indicated by star symbols. 

faster than the boson ground-state energy with the asymptotics ~ — l/r 2 J z . Therefore close 
to r s = the continuous transformation of bosons into fermions by varying of parameter v 
breaks down. This result is not surprising because for r s — >• we have no ideal boson gas. 

Altogether the curves in Figs. 1-3 demonstrate the capability of our analytical formula 
to describe the ground-state energy in a wide range of the density parameter r s and for 
the whole range of the anyon parameter v. The ground-state energies for Coulomb anyon 
systems would continuously change between the curves from Eq. ( 136|) and Eqs. ( |33i) and 
(I39p if v would change from to 1 (except the close region to r s — 0). 
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FIG. 2: The same as in Fig. 1, however, for range 0.2 <r s < 1.5. 

VI. CONCLUSIONS 

We have derived an approximate analytic formula for the ground-state energy of the 
Coulomb anyon gas. Starting from our previous results for the noninteracting 2D anyon gas 
confined in a harmonic potential, we have flattened out the confinement with simultaneous 
increasing of the particle number to obtain the thermodynamic limit at fixed density. We 
have generalized this result to describe the interacting 2D Coulomb anyon gas by introducing 
a function of the anyon parameters v and r s , that takes into account the Coulomb interac- 
tion and fractional statistics. We have determined this function by fitting to the analytic 
expressions for the ground-state energy of the classical electron crystal at very large r s , to 
that of the 2D Coulomb Bose gas at very small r s and to the HF energy (at high density) 
for spin-polarized 2D electrons. Our analytic formula applies to the full range of parameters 
v and r s and provides a convenient description of the ground-state of 2D Coulomb anyon 
systems. By comparison with HF and MC results from the literature we find significant 
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FIG. 3: Ground-state energies per particle vs anyon parameter v for various fixed values of r s . 
deviations for intermediate densities, which indicate shortcomings of these approaches. 
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